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Mathematical models for vaccination, waning 
immunity and immune system boosting: 
a general framework 
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Abstract 

When the body gets infected by a pathogen or receives a vaccine dose, the 
immune system develops pathogen-specific immunity. Induced immunity decays 
in time and years after recovery/vaccination the host might become susceptible 
again. Exposure to the pathogen in the environment boosts the immune system 
thus prolonging the duration of the protection. Such an interplay of within host 
and population level dynamics poses significant challenges in rigorous mathe¬ 
matical modeling of immuno-epidemiology. The aim of this paper is twofold. 
First, we provide an overview of existing models for waning of disease/vaccine- 
induced immunity and immune system boosting. Then a new modeling approach 
is proposed for SIRVS dynamics, monitoring the immune status of individuals 
and including both waning immunity and immune system boosting. We show 
that some previous models can be considered as special cases or approximations 
of our framework. 
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1 Introduction 

Models of SIRS type are a traditional topic in mathematical epidemiology. Classi¬ 
cal approaches present a population divided into susceptibles (S), infectives (I) and 
recovered (R), and consider interactions and transitions among these compartments 
|9]. Susceptibles are those hosts who either did not contract the disease in the past or 
lost immunity against the disease-causing pathogen. When a susceptible host gets in 
contact with an infective one, the pathogen can be transmitted from the infective to 
the susceptible and with a certain probability the susceptible host becomes infective 
himself. After pathogen clearance the infective host recovers and becomes immune for 
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some time, afterward he possibly becomes susceptible again (in certain cases one can 
talk of life-long immunity). The model can be extended by adding vaccination. Vac- 
cinees (V) are protected from infection for some time, usually shorter than naturally 
infected hosts. 

From the in-host point of view, immunity to a pathogen is the result of either active 
or passive immunization. The latter is a transient protection due to the transmission 
of antibodies from the mother to the fetus through the placenta. The newborn is thus 
immune for several months after birth ra- Active immunization is either induced by 
natural infection or can be achieved by vaccine administration [351 US- 

Let us first consider the case of natural infection. A susceptible host, also called 
naive host , has a very low level of specific immune cells for a pathogen (mostly a 
virus or a bacterium). The first response to a pathogen is nonspecific, as the innate 
immune system cannot recognize the physical structure of the pathogen. The innate 
immune response slows down the initial growth of the pathogen, while the adaptive 
(pathogen-specific) immune response is activated. Clonal expansion of specific immune 
cells (mostly antibodies or CTL cells) and pathogen clearance follow. The population 
of pathogen-specific immune cells is maintained for long time at a level that is much 
higher than in a naive host. These are the so-called memory cells and are activated 
in case of secondary infection (see Figure [H adapted from [6].). Memory cells rapidly 
activate the immune response and the host mostly shows mild or no symptoms [2]. 

Each exposure to the pathogen might have a boosting effect on the population of 
specific memory cells. Indeed, the immune system reacts to a new exposure as it did 
during primary infection, thus yielding an increased level of memory cells. Though 
persisting for long time after pathogen clearance, the memory cell population slowly 
decays and in the long run the host might lose his pathogen-specific immunity [37 . 

Vaccine-induced immunity works in a similar way as immunity induced by the 
natural infection. Agents contained in vaccines resemble, in a weaker form, the disease- 
causing pathogen and force a specific immune reaction without leading to the disease. 
If the vaccine is successful, the host is immunized for some time. Vaccinees experience 
immune system boosting and waning immunity, just as hosts recovered from natural 
infection do. In general, however, disease-induced immunity induces a much longer 
lasting protection than vaccine-induced immunity does [35]. 

Waning immunity might be one of the factors which cause, also in highly developed 
regions, recurrent outbreaks of infectious diseases such as measles, chickenpox and 
pertussis. On the other side, immune system boosting due to contact with infectives 
prolongs the protection duration. In a highly vaccinated population there are a lot of 
individuals with vaccine-induced immunity and few infection cases, as well as many 
individuals with low level of immunity. In other words, if a high portion of the popu¬ 
lation gets the vaccine, there are very few chances for exposure to the pathogen and 
consequently for immune system boosting in protected individuals. 

In order to understand the role played by waning immunity and immune system 
boosting in epidemic outbreaks, in the recent past several mathematical models were 
proposed. Few of these models describe only in-host processes during and after the 
infection [311 [16] . Many more models, formulated in terms of ordinary differential 
equations (ODEs), consider the problem only at population level, defining compart- 
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Figure 1: Level of pathogen-specific immune cells with respect to the time. The solid 
line represents the case of natural infection, the dotted line represents the immune 
status of a vaccinated host. Generation of memory cells takes a few weeks: once 
primary infection (respectively, vaccination) occurred, the adaptive immune system 
produces a high number of specific immune cells (clonal expansion). After pathogen 
clearance, specific immune cells (memory cells) are maintained for years at a level that 
is much higher than in a naive host. Memory cells are activated in case of secondary 
infection. 

ments for individuals with different levels of immunity and introducing transitions 
between these compartments mm- Vaccinated hosts or newborns with passive im¬ 
munity are often included in the model equations and waning of vaccine-induced or 
passive immunity are observed [32, 2D ID 2], 122, 3 2D]- 

To describe the sole waning immunity process, authors have sometimes chosen delay 
differential equation (DDE) models with constant or distributed delays |2T], [36l 8,0 38]. 
The delay represents the average duration of the disease-induced immunity. However, 
neither a constant nor a distributed delay allows for the description of immune system 
boosting. 

Models which include partial differential equations (PDEs) mostly describe an age- 
structured population [27, 20], 22] and consider pathogen transmission among the dif¬ 
ferent age groups (newborns, children, pupils, adults, ...). Rare examples suggest a 
physiologically structured approach with populations structured by the level of immu¬ 
nity, coupling within-host and between-hosts dynamics [23 E]- 

The goal of the present book chapter is twofold. On the one side, we found nec¬ 
essary to provide a comprehensive overview of previously published models for waning 
of disease/vaccine-induced immunity and immune system boosting (Sect. [2]). On 
the other side, in Sect. [2] we propose a new modeling framework for SIRVS dynamics, 
monitoring the immune status of individuals and including both waning immunity and 
immune system boosting. 
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2 Mathematical models for waning immunity and 
immune system boosting 

In the following we provide an overview on previous mathematical models for waning 
immunity and immune system boosting. We shall classify these models according to 
their mathematical structure (systems of ODEs, PDEs or DDEs). 

2.1 Systems of ODEs 

Mossong and coauthors were among the first to suggest the inclusion of individuals 
with waning immunity in classical SIRS systems [3Tj. Motivated by the observation 
that measles epidemics can occur even in highly vaccinated populations, the authors 
set up a model to study the waning of vaccine-induced immunity and failure of se¬ 
roconversion as possible causes for recurrent outbreaks. Their compartmental model 
includes hosts with the so-called “vaccine-modified measles infection” (VMMI) which 
can occur in people with some degree of passive immunity to the virus, including those 
previously vaccinated. Assuming that not all vaccinees are protected from developing 
VMMI, the authors classify vaccinees into three groups: immediately susceptible to 
VMMI (weak response), temporarily protected who become susceptible to VMMI due 
to waning of vaccine-induced immunity (intermediate response), and permanently pro¬ 
tected from VMMI (strong response). Infection occurs due to contact with infectious 
individuals (both regular measles infection and VMMI). The resulting compartmental 
model includes waning of vaccine-induced immunity but not of disease-induced immu¬ 
nity, nor immune system boosting. Similar to McLean and Blower [28], Mossong et al. 
define a parameter 0 to describe the impact of the vaccine: if 0 < 1, then vaccine 
failure is possible. Analytical results in [31] show that the main effect of VMMI is to 
increase the overall reproduction number of the infection. 

Inspired by Mossong’s work, in 2003-2004 Glass, Grenfell and coauthors mmm 
proposed modifications and extensions of the system in [31]. The basic model is similar 
to the ODE system in [31], with a group of subclinical cases which carry the pathogen 
without showing symptoms m In addition, the distribution of antibody levels in 
immune hosts (included in the ODEs coefficients) and immune system boosting are 
introduced: the average antibody level in an immune host increases due to contact 
with infective or subclinical hosts. This model was used to fit measles data in Eng¬ 
land [123 • I n [13] the basic model was extended to consider measles transmission in a 
meta-population with N patches. 

Immune system boosting in vaccinees was further studied in [3]. In this work two 
models are introduced. In the first one vaccinees are separated from non-vaccinated 
hosts. Both groups of individuals are classified into susceptible, infective and immune, 
but in contrast to the models in [141 T3I [121 3T, there is no compartment for subclinical 
cases. Non-vaccinated hosts do not undergo immune system boosting. For vaccinated 
hosts the authors include a so called “self-boosting” of vaccine, so that contact with 
infectives moves susceptible vaccinees to the immune vaccinated compartment. The 
second model extends the first one with a new compartment for hosts with waning 
immunity (W). These can receive immune system boosting due to contact with infec¬ 
tives or move back to the susceptible compartment due to immunity loss. Numerical 
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simulations show possible sustained oscillations. The SIRWS system was partially an¬ 
alyzed by Dafilis et ah 

Heffernan and Keeling |16j proposed an in-host model to understand the behavior 
of the immune system during and after an infection. Activation of immune system ef¬ 
fectors and production of memory cells depend on the virus load. When not stimulated 
by the virus, the number of activated cells decays (waning immunity). Vaccination is 
simulated by changing the initial conditions for the virus load. Numerical simulations 
show that the number of infected immune system cells in a vaccinated patient reaches 
approximately half of what is reached in a patient who undergoes natural infection. 
In turn, the level of immunity gained after one dose of vaccine is the same as the level 
observed in a measles patient 4 years after natural infection. The in-host model in 
[16] was extended by the same authors to a population model (SEIRS) with waning 
immunity and immune system boosting m ■ In contrast to classical SEIRS models, 
the class R refers here to individuals protected by short-term immune memory, while 
the class S refers to those individuals who have lost this short-term protection and may 
experience immune system boosting. Each compartment is classified according to the 
level of immunity, which can be related to the number of memory cells. Newborns 
are recruited into the susceptible class So (lowest level of immunity). During exposure 
and infection the host does not change his level of immunity, that is, transition occurs 
from Sj to Ej to Ij for each j E N. Hosts in S and R experience waning immunity 
and transit from Sj to Sj_ i (respectively from Rj to Rj-i). Immune system boosting 
is due to recovery from infection and is incorporated into the equations with transi¬ 
tion terms from Ij to Rk, with k > j. The resulting large system of ODEs, with a 
very high number of parameters, is quite hard to approach from an analytical point of 
view, hence the authors make use of numerical simulations to investigate the long term 
behavior. A somehow simplified version of the ODE system in di was proposed by 
Reluga et al. in [32], A similar large system of ODEs was introduced by Lavine et al. 
ra, extending the SIRWS model in mm, by including several levels of immunity 
for immune hosts (R) and hosts with waning immunity (W), as well as age classes 
for all compartments. The authors claim that the model can explain several observed 
features of pertussis in US, in particular a shift in the age-specific incidence and the 
re-emergence of the disease in a highly vaccinated population. 

2.2 System of DDEs 

Delay models with constant or distributed delay have been introduced to describe 
waning of disease-induced or vaccine-induced immunity. A simple SIRS system with 
constant delay is given by 

S(t) = Ml - S(t)) - 4>S(t)f(m) + 7 I(t - T)e-“ r 

i(t) = ts(t)f(m)-(p+7)m (i) 

R(t) = 7 1{t ) — /ui?(t) —7 lit — r)e _/iT . 

This model was studied by Kyrychko and Blyuss ra. who provided results on exis¬ 
tence, uniqueness and non-negativity of solutions, linear and global stability of the 
disease-free equilibrium, as well as global stability of the unique endemic equilibrium. 
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A special case of ([I]) was considered some years later by Taylor and Carr |36j . An 
extension of system ([I]) with distributed delay was proposed in [8] and shortly after in 

0 - 

A more general model with distributed delay and vaccination was proposed by 
Arino et al. in [3], Their system includes three compartments (susceptible, infective 
and vaccinated hosts) in a population which remains constant in time. Vaccine-induced 
immunity might be only partial, resulting in vaccinated individuals becoming infec¬ 
tive. Systems of ODEs or DDEs can be obtained from the general model by a proper 
choice of the kernel (see also [T9l [I8] h 

Recently, Yuan and Belair proposed a SEIRS model with integro-differential equa¬ 
tions which resembles the systems in [MH]. The probability that an individual stays 
in the exposed class (E) for t units of time is P(t), hence, 

Similarly, Q{t) is the probability that an individual is immune t units of time after 
recovery, thus 

R{t) = [ 7 1{u)e- b{t ~ u) Q{t - u) du. 

Jo 

For a certain choice of the probabilities P and Q, the problem can be reduced to a 
system with one or two constant delays. The authors show existence of an endemic 
equilibrium and boundedness of solutions in a positive simplex. For the system with 
one constant delay, results for existence of a global attractor as well as the proof of 
persistence of the disease in case Ro > 1 are provided. 

2.3 Systems of PDEs 

Structured populations in the context of waning immunity and immune system boost¬ 
ing have been motivated in different ways. Often the structure can be found in the 
biological age (2ZIG2H EH M\, and is used to observe disease transmission among ba¬ 
bies, children, adults and seniors. Only few works suggest models for physiologically 
structured populations mm- 

McLean and Anderson mm proposed a model for measles transmission which in¬ 
cludes a compartment for babies protected by maternal antibodies. Indeed, mothers 
who have had measles or have been vaccinated transfer measles immunity to the baby 
through the placenta. For several months after birth (ca. 2 months if the mother 
was vaccinated, ca. 4 months if she had the disease |2B]) the baby is still protected 
by maternal antibodies and should not be vaccinated. The model by McLean and 
Anderson |23 considers only waning of maternally induced immunity in the context of 
measles infection. Few years before McLean, Katzmann and Dietz m proposed a bit 
more general model, which includes also waning of vaccine-induced immunity. In both 
cases, the age structure was used to determine the optimal age for vaccination. A com¬ 
partment for adult hosts with waning immunity who can also receive immune system 
boosting was introduced only years later by Rouderfer et al. [33]. A further determin¬ 
istic system of ODEs for maternally induced immunity in measles was proposed in [29]. 
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Different is the approach when physiologically structured populations are considered. 
Martcheva and Pilyugin [25] suggest an SIRS model in which infective and recovered 
hosts are structured by their immune status. In infective hosts the immune status 
increases over the course of infection, while in recovered hosts the immune status de¬ 
cays at some non-constant rate. When the immune status has reached a critical level, 
recovered hosts transit from the immune to the susceptible compartment. 

A general framework for SIRS systems, modeling waning immunity and immune 
system boosting, and combining the in-host perspective with the population dynamics, 
was proposed in [ 6 ]. 

3 A general modeling framework 

In this section we extend the model in [ 6 ] to include vaccine-induced immunity. As in 
[25116], we couple the in-host with the between-hosts dynamics, focusing on the effects 
of waning immunity and immune system boosting on the population dynamics. In 
contrast to the models proposed in [T7] E2], we shall maintain the number of equations 
as low as possible. The resulting model (VI) is a system of ODEs coupled with two 
PDEs. The ODE systems in [31] QT] [3] [5] |30], as well as extensions of the DDEs 
systems in [ 361 l38 |, can be recovered from our modeling framework. 

Setting up our model we do not restrict ourselves to a particular pathogen. The 
model (VI) can be adapted to several epidemic outbreaks (e.g. measles, chickenpox, 
rubella, pertussis) by ad-hoc estimating coefficients from available experimental data 

[211 HI E3!- 

3.1 Model ingredients 

3.1.1 Originally susceptible and infectives hosts 

Let S(t ) denote the total population of originally susceptible hosts. These are sus¬ 
ceptible individuals which have neither received vaccination nor have been infected 
before. Newborns enter the susceptible population at rate b(N), dependent on the to¬ 
tal population size N. For simplicity we assume that the natural death rate d > 0 does 
not depend on N. Assume that b : [0, oo) —>■ [0,6+], N (-)■ b(N), with 0 < 6 + < oo, 
is a nonnegative function, with 6(0) = 0. Finally, assume that in absence of disease- 
induced death there exists an equilibrium N* such that b(N*) = dN*. 

Let I{t ) denote the total infective population at time t. Infection of susceptible 
individuals occurs by contact, at rate /3I/N. Infected hosts recover at rate 7 > 0. 
When we include disease-induced death at rate dj > 0, the equilibrium N* satisfies 

b(N*) = dN* + d T I*. 


3.1.2 Immune individuals 

Let us denote by r(t, z) the density of recovered individuals with disease-induced im¬ 
munity level z G [z rnin , z max \ at time t. The total population of recovered hosts is given 
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by 

f*Zmax 

I r(t , z) dz. 

Zmin 

The parameter z describes the immune status and can be related to the number of 
specific immune cells of the host. The value z max corresponds to maximal immunity, 
whereas z min corresponds to low level of immunity. Individuals who recover at time 
t enter the immune compartment with maximal level of immunity z max ■ The level 
of immunity tends to decay in time and when it reaches the minimal value z m i n , the 
host becomes susceptible again. However, exposure to the pathogen can boost the 
immune system from z £ [z m i n , z max ] to any higher status. It is not straightforward 
to determine how this kind of immune system boosting works, as no experimental 
data are available. Nevertheless, laboratory analysis on vaccines tested on animals 
or humans suggest that the boosting efficacy might depend on several factors, among 
which the current immune status of the recovered host and the amount of pathogen 
he receives mm- Possibly, exposure to the pathogen can restore the maximal level 
of immunity, just as natural infection does |JIj. 

Let p(z, z), z > z, z, z £ M denote the probability that an individual with immunity 
level 5 moves to immunity level z , when exposed to the pathogen. Due to the definition 
of p(z , z), we have p(z, z) £ [ 0 ,1], z > z and 

p(z, z) = 0 , for all z < z. 

As we effectively consider only immunity levels in the interval [z min , z max ], we set 
p(z, z) = 0, for all 5 £ (-oo, z min ) U ( z max , oo). 



Then we have 

/ OO rZmax 

p(z,z)dz = / p(z, z)dz = 1 , for all 5 £ [z min , z max \. 

-oo J z 

Exposure to the pathogen might restore exactly the immunity level induced by the 
disease ( z max ). In order to capture this particular aspect of immune system boosting, 
we write the probability p(z, z) as the combination of a continuous (p 0 ) and atomic 
measures (Dirac delta): 


p{z, z) = c max (z)S(z max ~z) + c 0 {z)po(z, z) + ci{z)S(z - z), 


where 


• C max ■ [ Zmin, z m ax\ [0,1], y c max (y), is a continuously differentiable function 
and describes the probability that, due to contact with infectives, a host with 
immunity level y boosts to the maximal level of immunity z max . 

• c 0 : [zmin, z m ax] [ 0 , 1 ], y c 0 (y), is a continuously differentiable function and 
describes the probability that, due to contact with infectives, a host with immu¬ 
nity level y boosts to any other level z £ (y,z max ), according to the continuous 
probability po(z, y ). 

• ci(y) = 1 — c max {y ) — co(y) describes the probability that getting in contact with 
infectives, the host with immunity level y £ [z m in,z max ] does not experience 
immune system boosting. 

8 


The immunity level decays in time at some rate g(z) which is the same for all 
recovered individuals with immunity level z. In other words, the immunity level z 
follows 


d 

dt 


z{t) = 9(z), 


with g : [z mm , z max \ —> (0, K g ], K g < oo continuously differentiable. The positivity 
of g(z) is required from the biological motivation. Indeed, if g(z) = 0 for some value 
5 G [ z min , z max \, there would be no change of the immunity level at z, contradicting the 
hypothesis of natural decay of immune status. In absence of immune system boosting, 
we have that 



1 

9{x) 


dx 


is the time a recovered host remains immune (see El)- 


3.1.3 Vaccination 


We structure the vaccinated population by the level of immunity as well. Let v(t,z ) 
be the density of vaccinees with immunity level z G [z mini z max ] at time t. The total 
population of vaccinated hosts is given by 


V(t) = 


v(t , z) dz. 


Vaccination infers a level of immunity z vax , which is lower than the level of immunity 
after natural infection: z max > z vax > z min [35]. As in recovered individuals, the 
level of immunity of a vaccinated host tends to decay in time and when it reaches the 
minimal value z m i n , the host becomes susceptible again. However, also in vaccinated 
hosts, exposure to the pathogen can boost the immunity level z G [z m i n , z vax ] to any 
higher value in [z mm , z max \. Immune system boosting is described by the probability 
p(z, z), as in recovered hosts. We consider the possibility that exposure to the pathogen 
boosts the immune system of a vaccinated individual to z G (z vax , z max J. Vaccinated 
hosts with ^ G (z vax ,z max \ have an immune status which can be compared to the one 
of hosts who recovered from natural infection. 

It is reasonable to assume that in vaccinated individuals the immunity level decays 
in time at the same rate g , as in hosts who underwent natural infection. In absence 
of exposure to the pathogen (hence in absence of immune system boosting), the time 
that a vaccinee remains immune is shorter than the time a recovered host does: 

PZvax f*Zmax 

Jz min g(x) dx < Jz min g(x) dx ' 


Let us define the vaccination rate at birth a > 0. We assume that originally susceptible 
(adult) individuals get vaccinated at rate 4> > 0. 


3.1.4 Becoming susceptible again 

In absence of immune system boosting both disease-induced and vaccine-induced im¬ 
munity fade away. Individuals who lose immunity either after recovery from infection 
or after vaccination, enter the class S 2 of susceptible individuals who shall not get a 
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new dose of vaccine. A host who had the disease or got vaccination relies indeed on 
the induced-immunity and is not aware of the fact that his level of immunity might 
have dropped below the critical immunity threshold. 


We denote by 5' 2 (t) the population at time t of susceptible hosts who are not go¬ 
ing to receive vaccination. 

3.2 Model equations 

In view of all what we have mentioned above, we can easily write down the equations 
for the compartments S, I and S'2. Let initial values S'(O) = S° > 0, 1(0) = 1° > 0 
and 52(0) = S'® > 0 be given. The population of originally susceptible individuals is 
governed by 


S(t) = b{N(t))(l-a)~ (j>S{t) ~ dS(t), 


birth 


vaccination 


infection 


( 2 ) 


death 


whereas hosts who become susceptible due to immunity loss follow 

Mt)m 


S 2 (t) = -P 


infection 


dS2(t) + Ar + 


A 


v. 


death. immunity loss immunity loss 
after recovery after vaccination 


The term A# (respectively Ay), which represents transitions from the immune (re¬ 
spectively, the vaccinated) compartment to the susceptible one, will be specified below 
together with the dynamics of the recovered (respectively, vaccinated) population. 


Both kinds of susceptible hosts can become infective due to contact with infective 
hosts: 


■, . S(t)I(t ) ) . . . . . . 

{t)= + 13 N(t ) ~ 2 ^ -dm- dim • 

^ v ^ ^ v ^ recovery natural disease-induced 

infection of s infection of s 2 death death 


(3) 


To obtain an equation for the recovered individuals, structured by their levels of im¬ 
munity, one can proceed similarly to size structured models or as it was done for the 
immune population in [6]. The result is the following PDE. Let a nonnegative initial 
distribution r(0,z) = ip(z), z G [z min , z max ] be given. For t > 0,z G [z min ,z max \ we 
have 


^r(t, ~ z)) = -dr(t, z) + 


p(z, x)r(t, x) dx 


r(t,z)P 


m 

N(tY 


(4) 


with the boundary condition 


Q^Zmax )r(t, z ma x) qd(^) T P 


m 

N(t ) 


P(Zmax,x)r(t,x) dx. 


(5) 
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Equation 03]) expresses the rate of change in the density of recovered individuals ac¬ 
cording to immune level due to natural waning, mortality, and boosting. The boundary 
condition OS]) includes newly recovered individuals as well as those recovered individ¬ 
uals, who just received a boost which elevated their immune system to maximal level. 

Next we shall consider the vaccinated population. Again, by structuring this group 
according to immunity level, one has the PDE 


J^(t, Z ) = ffZ ( 9(z)v(t , z)) - dv(t, Z ) + /3-^y 

m 


p(z, x)v(t, x) dx 


( 6 ) 


- v(t, y + S(z - z vax ) ( 4>s(t ) + ab(N(t ))), 


and 


9[z , max)'V(ti Z rnax ) ft 


m 

m 


p(z max ,x)v(t,x) dx, 


(7) 


provided with a nonnegative initial distribution v(Q,z) = ipv{z), z G [z m i n , z max \. Ob¬ 
serve that newly vaccinated hosts do not enter the vaccinated population at z max , but 
at the lower value z vax , which is expressed in equation OS]) as an impulse at z — z vax 
by the term with the Dirac delta S(z — z vax ) ■ 

ft becomes evident that the quantity Ar, initially introduced in the S2 equation 
to represent the number of hosts who experienced immunity loss, is given by the 
number g{z mi , n )r{t, z rnin ) of immune hosts who reached the minimal level of immunity 
after recovery from natural infection. Similarly, Ay is the number g{z mi , n )v{t, z min ) of 
vaccinated hosts who reached the minimal level of immunity. Hence we have 


s 2 (t) 


swm 
m , 

infection 


^£ 2(^0 + T g{^ min)v(di ^min ) • 

death A v 


( 8 ) 


In the following we refer to the complete system 02]) - OH) as to model (VI). 


4 Connection to other mathematical models 

4.1 Connection to ODE models 

As it was shown in |6j for a simpler problem, model (VI) can be reduced to a system 
of ODEs analogous to those proposed in mmmmm [30 |. The connection be¬ 
tween model (VI) and the ODE system is given by the method of lines , a technique in 
which all but one dimensions are discretized [33]. In our case, we shall discretize the 
immunity level (z) and obtain a system of ODEs in the time variable. 

Let us define a sequence with hj := z J+ \ — Zj > 0, for all j 6 N. To 

keep the demonstration as simple as possible, we choose a grid with only a few points, 
Z\ := z min < Zw z vax < Zp < z max and for simplicity (or possibly after a rescal¬ 
ing) assume that hj = 1 for all j. We define the following subclasses of the im¬ 
mune/vaccinated population: 
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R F (t) := r(t,Zp)i immune hosts with high level of immunity at time t. As their 
immunity level is quite high, these individuals do not experience immune system 
boosting. Immunity level decays at rate // := g(z F ) > 0. 


• Rw(t) f(t,Zw), immune hosts with intermediate level of immunity at time t. 
These individuals can get immune system boosting and move to Rf■ Immunity 
level decays at rate v := g(zw) > 0 . 

• Rc(t) := r(t,z m i n ), immune hosts with critically low level of immunity at time 
t. With probability 8 boosting moves Rc individuals to R,\y (respectively, with 
probability (1 — 8) to Rf)- Immunity level decays at rate a := g(z min ) > 0 . 
If they do not get immune system boosting, these hosts move to the class S 2 
(become susceptible again). 

• Vr(£) := v(t,z F ), vaccinated hosts who thanks to immune system boosting 
gained a very high level of immunity at time t. These individuals do not ex¬ 
perience immune system boosting. Immunity level decays at rate //. 

• Vo (t) := v(t, zw), vaccinated individuals at time t with maximal vaccine-induced 
immunity. This class includes new vaccinees. If their immune system gets 
boosted hosts move to Vr. Immunity level decays at rate v. 

• Vc{t) := v(t,z m i n ), vaccinees with critically low level of immunity at time t. 
With probability £ boosting moves Vc hosts to Vo and with probability (1 — £) 
to Vr. Immunity level decays at rate cr. If they do not receive immune system 
boosting, Vc hosts move to S 2 - 

To show how the PDE system can be reduced to a system of ODEs by means of the 
method of lines, we consider a simple example. Let us neglect immune system boosting 
for a moment. Then the PDE for r(t,z ) in model (VI) becomes 


— r(t, z) = —(g(z)r(t, z)) - dr(t, z), z e [z min , z max ], 


(9) 


with boundary condition R Zmax (t) := r(t,z max ) = 7 1(t)/g(z max ). Using forward ap¬ 
proximation for the ^-derivative in (JUJ), we obtain, e.g., for Rp(t ) the following differ¬ 
ential equation: 



g(z max )r(t, z max ) - g{z F )r(t, z F ) _ ^ 


max 


max 



9( Z max)Rzmax{t) gR F (t) dR F (t ) 
= 7 1(t) — {n + d)R F (t). 
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Analogously one can find equations for R w , Rq, V r , V 0 and Vc- Altogether we obtain 
a system of ordinary differential equations in which a linear chain of ODEs replaces 
the PDEs for the immune and the vaccinated class: 


S(i) = (1 - a)b(N(t)) - m) - /Hvyr " dS{t) 

m = + S M) - n+ d + d,)i(t) 

R F (t) = 7 /(f) — jxR F (t) — dR F (t) 

Rw(t) = I^Rf(P) — vRw{t) — dR w (t) 

Rc(t) = uR w (t) — aRc{t ) — dR c (t ) 

V R (t) = ~nV R {t) - dV R {t ) 

V Q {t) = 4bS{t ) + ab(N(t)) + pV R {t) - uV 0 (t) - dV 0 (t ) 

Veit) = vV 0 (t) - crV c (t ) - dV c (t) 

7 (t) = - dS 2 (t) + a(R c (t) + V c (t )). 

The method of lines can be applied to the full model (VI) as well [6]. To this purpose 
it is necessary to discretize the boosting probability p(z, z) (this is expressed by the 
parameters £ and 9 below). Incorporating the boosting effect, the result is the following 
system of ODEs. 


S(t) 

i(t ) 

R F (t) 

Rw(t) 

Rc(t ) 
Vn(t) 

Vo (t) 
Veit) 

S 2 (t) 


(1 - a)b(N(t)) - 4>S(t) - - dS(t) 


N(t) 


m 


PjjTk(S(t) + S 2 (t))-('y + d + d I )I(t) 
ll{t ) - nR F {t) - dR F {t) + pj^ (. R w (t ) + (1 - 0)Rc(t)) 
pR F (t) — uR w (t) — dR w (t ) + P (0Rc(t) ~ Rw(t)) 

vR w (t) - aR c (t ) - dR c (t ) - P^Rc(t) 

^N§) {Vo{t) + (1 “ ®Mt)) - pV R (t) - dv R (t) 

pS{t) + ab{N{t)) + pV R (t) - uV 0 {t) - dV 0 (t) + pj^L ^ Vc (t) - V Q (t)) 


= vV 0 (t) - aV c (t) - dV c (t) - P 


V c (t)I(t) 
Nit) 


' P ^N “ dSM + a{Rc(t) + VcW) ' 


The linear chain of ODEs provides a rough approximation of the PDEs in model (VI). 
Indeed, with the method of lines we approximate the PDE dynamics considering only 
changes at the grid points ( z m i n , Zw, z F ), whereas the dynamics remains unchanged 
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in each immunity interval [_Zj, z 3+ \ ]. We consider as representative point of the interval 
the lowest boundary z 3 - for this reason we do not have a differential equation for 
R , (t) or V z (t). 

z max \ / z max \ / 


4.2 Connection to DDE models 

Delay models with constant delay can be recovered from special cases of model (VI). 
We show here how to obtain the classical SIRS model with delay studied in [36], or 
extensions thereof. 


In the following we neglect boosting effects and vaccination. Further we do not distin¬ 
guish between originally susceptibles and host who have lost immunity, hence w.r.t. 
model (VI) we identify the classes S and £> 2 . From our assumptions, the disease- 
induced immunity lasts for a fix time, r > 0 years, given by 



1 

g{%) 


dx = T. 


We can express the total immune population at time t as the number of individuals 
who recovered in the time interval [t — r, t\, 

R(t) = 7 f I{y)e~ d{t - y) dy = 'y [ I(t- x)e~ dx dx. 

Jt-T JO 


Differentiation with respect to t yields 

R(t) = 7 /(t) — 7 1(t — r)e~ dT — dR(t). (10) 


On the other side, we have the definition in terms of distribution of immune individuals, 



Differentiate the last relation and compare with (fTOl) : 

9 (Zmax)^ (tj ^max ) 0 / -^(^')> 5 , (^ / min)^’(0 ^min ) v)e 


This means that individuals with maximal level of immunity are those who recover 
from infection. If a host who recovers at time t\ survives up to time t± + r, he exits 
the R class and enter S. In turn, we find a delay term in the equation for S too, and 
have a classical SIRS model with constant delay 

S(t) = b(N(t)) - - dS(t) + 7 1(t - r)e~ dT 

m = p^§^--(i+ d + d i)m 

R(t) = 7 1(t) —7 1(t — r)e~ dT — dR(t ), 
which was studied by Taylor and Carr [36]. 
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Now we can include again vaccination and the class S 2 as in the general model (VI). 
We assume that vaccine-induced immunity lasts for a time r v > 0, 


To, := 


g(x] 


dx < 


g(x] 


■ dx =: t. 


With similar arguments as for the immune population, we obtain the relations 


g^vax )v(t 1 Zvax ) = ab(N(t )) + 4>S(t), 
g{z m in)v(t, z min ) = ( ab(N(t - Ty)) + 4>S(t- Ty)) e ~ dTv , 

and find a system with two constant delays 

S(t) = (1 - a)b(NX) - <t>s(t) - _ dS (t) 

U) = ' /^ m + ^2(0) - (7 + d + 

R(t) = 7 1 [t) — 7 1(t — r)e~ dT — dR(t) 

V (f) = ab(N(t)) + (/>S(t) — (ab(N(t — t v )) + (f>S(t — t v )) e~ dTv — dV (f) 
s 2 (t) = - dS 2 (t) + 7 1(t - r)e-*- 

+ ( ab(N(t - T v )) + 4>S(t - Ty)) e~ dTv . 
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